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ABSTRACT 



The rate of mass accumulation due to galaxy merging depends on the mass, 
density, and velocity distribution of galaxies in the near neighborhood of a host 
galaxy. The fractional luminosity in kinematic pairs combines all of these ef- 
fects in a single estimator which is relatively insensitive to population evolu- 
tion. Here we use a k-corrected and evolution compensated volume-limited 
sample having an R-band absolute magnitude of M^' e < —19.8 + 51og/i mag 
drawing about 300 redshifts from CFGRS and 3000 from CNOC2 to measure 
the rate and redshift evolution of merging. The combined sample has an ap- 
proximately constant co-moving number and luminosity density from redshift 
0.1 to 1.1 (Vt M = 0.2, Qa = 0.8); hence, any merger evolution will be domi- 
nated by correlation and velocity evolution, not density evolution. We identify 
kinematic pairs with projected separations less than either 50 or 100 h~ l kpc 
and rest-frame velocity differences of less than lOOOkms -1 . The fractional lu- 
minosity in pairs is modeled as /l(Ai> , r p , M^ e )(l + z) mL where are 
[0.14±0.07,0±1.4] and [0.37 ±0.7, 0.1 ±0.5] for r p < 50 and 100/i _1 kpc, respec- 
tively (Q M — 0.2, f2 A = 0.8). The value of m L is about 0.6 larger if A = 0. To 
convert these redshift space statistics to a merger rate we use the data to derive a 
conversion factor to physical space pair density, a merger probability and a mean 
in-spiral time. The resulting mass accretion rate per galaxy (Mi, M 2 > 0.2M*) is 
0.02 ± 0.01(1 ±2) 01±a5 M* Gyr -1 . Present day high-luminosity galaxies therefore 
have accreted approximately 0.15M* of their mass over the approximately 7 Gyr 
to redshift one. Since merging is likely only weakly dependent on host mass, the 
fractional effect, 5M/M ~ 0.15M*/M, is dramatic for lower mass galaxies but is, 
on the average, effectively perturbative for galaxies above M*. 
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Subject headings: cosmology: large scale structure, galaxies: evolution 

1. Introduction 

Merging is a fundamental mode of stellar mass addition to galaxies. Moreover, 
merging brings in new gas and creates gravitational disturbances that enhance star 
formation or fuel a nuclear black hole. The general process of substructure infall may be 
the rate fixing process for the buildup of a galaxy's stars and consequently may largely 
regulate its luminosity history. Gravitational forces on relatively large scales dominate 
merger dynamics which allows direct observation of the mechanism, although with the 
considerable complication that dark matter dominates the mass. N-body simulations 
(Toomre & Toomre 1972; Barnes & Hernquist 1992) give the detailed orbital evolution, 
morphological disturbances and eventual outcomes of the encounters of pairs of galaxies. 

The purpose of this paper is to estimate the rate of mass gain per galaxy due to 
mergers over the redshift zero to one interval. Our primary statistic is the fractional 
luminosity in close kinematic pairs, which is readily related to n-body simulations and 
sidesteps morphological interpretation. This approach provides a clear sample definition 
which is closely connected to the large scale dynamics of merging. In common with all 
merger estimates it requires an estimate of the fraction of the pairs that will merge and a 
mean time to merger. 

The number of kinematic pairs is proportional to the volume integral at small 
scales of the product of two-point correlation function, £, and the luminosity function 
(LF). The high luminosity galaxies appear to be evolving purely in luminosity 
(Lilly et al. 1995; Lin et al. 1999), which can be easily compensated. The measured 
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evolution of £ suggests that the density of physical pairs should not vary much with redshift, 
(1 + z) 0±1 (LeFevre et al. 1996; Carlberg et al. 1997; Carlberg et al. 1999). This inference 
is in notable contrast with the pair counts or morphological typing approaches to merger 
estimation (Zepf & Koo 1989; Carlberg Pritchet & Infante 1994; Yee & Ellingson 1995; 
Patton et al. 1997; LeFevre et al. 1999), which suggest that merging rate by number varies 
as (1 + z) 3±1 . HST photometric pairs, with no redshift information leads to a dependence 
of (1 + Z ) l - 2±QA (Neuschaefer et al. 1997). 

In the next section we combine the Caltech Faint Galaxy Redshift Survey (CFGRS) 
and the Canadian Network for Observational Cosmology field galaxy survey (CNOC2) 
from which we construct evolution compensated, volume-limited, subsamples. In Section 
3 we measure the fractional luminosity in 50 and 100/iT 1 kpc companions as a function of 
redshift. The CNOC2 sample is used in Section 4 to relate this wide pair sample to a close 
pair sample which is more securely converted into a mass merger rate. Section 5 discusses 
our conclusions. We use H = 100/ikms -1 Mpc -1 , VL M = 0.2 in open and flat cosmologies. 

2. The CFGRS and CNOC2 Volume-Limited Samples 

The CFGRS sample of the HDF plus flanking fields is discussed in detail elsewhere 
(Hogg et al 1999; Cohen et al. 2000). We use the high coverage subsample lying within 
a 240 arcsecond radius circle, with a center located at \2 h 36 m 50 s and 62° 12' 55" 
(J2000). The computed magnitude selection function, s(m R ), (in Cousins R) is accurately 
approximated as a constant 90% spectroscopic completeness for m R < 22.8 mag with 
a linear decline to 19% at m R < 23.4 mag, our sample limit. The magnitude weight is 
l/s(m R ). The CFGRS k-corrections and evolution compensation are here approximated as 
k(z) = Kz mag from the tables of Poggianti (1997). For galaxies that Cohen et al. (2000) 
classify as "E" (emission), K = 1.0, "A" types have K = 2.0 and all types have K = 1.7 
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The CN0C2 selection weights and k-corrections are discussed in Yee et al. (2000). . 
The evolution of the luminosity function is approximated as a uniform M* (z) = M* — Qz, 
with (5 — 1, (Lin et al. 1999) which we use over the entire CNOC2-CFGRS redshift range. 

The kinematic pair fraction is directly proportional to the mean density of the sample 
and is therefore sensitively dependent on correcting to a complete and uniform sample 
(Patton et al. 2000a). The most straightforward approach is to impose a strict volume 
limit. For our primary sample we will limit the CFGRS and the CNOC2 samples at 
Mft £ = —19.8 + 51og/i mag, which yields volume-limited samples of about 300 CFGRS 
galaxies between redshift 0.3 to 1.1 and 3000 CNOC2 galaxies between 0.1 to 0.5. The 
volume density of the sample is approximately constant at 1.2 x 10~ 2 /i 3 Mpc -3 over the 
entire redshift range for il\ = 0.8 but rises roughly as (1 + z) os for il\ = 0. Both the 
CFGRS and the CNOC2 surveys are multiply masked, which minimizes the effects of slit 
crowding, however there is still a measurable pair selection function. The CNOC2 catalogue 
has about a 20% deficiency of close angular pairs. We model the measured angular pair 
selection weight as, w{6) = [1 + a s tanh(6 l /^ s )] _1 , where [a s ,0 s ] is [0.5,5"] for the CFGRS 
sample and [—0.3, 10"] for the CNOC2 sample with typical pair corrections being 10%. 

3. The Pair Fractional Luminosity Fraction 

The preferred choice of pair statistic depends on the application (Patton et al. 2000a). 
Here we are primarily interested in the impact of merging on galaxy mass increase, for 
which the k-corrected, evolution compensated R luminosity is a stand-in. The rate of 
merging per galaxy depends on the density of galaxies in the near neighborhood and their 
velocity distribution. As a practical redshift space estimator, we compute the fractional 
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luminosity in close kinematic pairs, 

f f*IA„mi« .max M k,es Ej Ei#,< A^",< rr WjW^jj^ 

f L {z\Av ,r ,M R ) = — , (1) 

where the weights, Wi, allow for the magnitude selection function. Note that the ij and 
ji pairs are both counted. The ratio has the benefit of being fairly stable for different 
luminosity limits, self-normalizing for luminosity evolution, identical to a mass ratio for a 
fixed M/L population. For an unperturbed pair luminosity function it is mathematically 
identical to the N c of Patton et at. (2000a) although constructed out of somewhat different 
quantities. The two parameters, Au max and r™ ax , are chosen on the basis of merger 
dynamics and the characteristics of the sample. The rate of mass increase per galaxy is 
calculated from this statistical estimator using a knowledge of merger dynamics and the 
measured correlations and kinematics of galaxy pairs in the sample. 

The mean fractional pair luminosity, based on 18 CFGRS pairs and 91 CNOC2 
pairs, with Av < 1000 km s" 1 and 5 < r p < 50/i _1 kpc pairs is displayed in Figure 1. 
These kinematic separation parameters are larger than is suitable for reliably identifying 
"soon-to-merge" pairs. However, they provide a statistically robust connection to those pairs 
and take into account the lower velocity precision and sample size of the CFGRS relative 
to CNOC2. The errors are computed from the pair counts, n p 1 ^ . The measurements of 
f L (z) in Figure 1 are fit to f L (Av, r p , M r fee )(l + z) mL , finding [f L ,m L ] of [0.14 ± 0.07, ± 1.4] 
for r p < 50/r 1 kpc pairs and [0.37 ± 0.7,0.1 ± 0.5] for r p < lOO/i^kpc, both for 
Vl M = 0.2, fi A = 0.8. The increase with r p of f L is consistent with a 7 = 1.8 two-point 
correlation function. If VLm = 0.3, = 0.7, then at 100/i -1 kpc rises by about 0.05 
whereas if Vl M = 0.2, = 0.0, then m L = 0.50. The increase is largely as a result of the 
rise in the implied co-moving sample density over this redshift range. 

The merger probability of a kinematic pair depends sensitively on the pairwise velocity 
dispersion, <t 12 , of galaxies. The model pairwise velocity distribution is computed as the 
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convolution of the correlation function with the distribution of random velocities. The infall 
velocities are negligible at these small separations and we will assume that the peculiar 
velocities are drawn from a Gaussian distribution. The measured fraction of the CNOC2 
pair sample with velocities smaller than some Av, normalized to the value at 1000 km s -1 , 
is displayed in Figure 2. The 50/i _1 kpc wide pairs limited at —19.5 mag are plotted as open 
squares, the 20/i -1 kpc pairs limited at —18.5 and —19.5 mag are plotted as octagons and 
diamonds, respectively. The upper curve assumes that a±2 is 200 kms -1 and the lower one 
300 kms -1 , which approximately span the data. 

4. Merger Rate Estimation 

The merger rate is best estimated from very close kinematic pairs, 20h~ 1 kpc or less, 
about half of which are physically close and have significant morphological disturbance 
(Patton et al. 2000a). However the fraction of galaxies in such close pairs is about one 
percent, giving poor statistics. Since the number of pairs increases smoothly as r 3 ~ 7 , where 
7 is the slope of the small scale correlation function, we can use pairs at somewhat larger 
separations as statistically representative of the close pairs, however we prefer to stay within 
the radius of virialized material around a galaxy over our redshift range, which is no larger 
than about lOO/i -1 kpc. The mass accretion rate from major mergers is therefore estimated 
as, 

Km = ^f L (Av,r p ,z)C zs (^v,l)F(v < v mg )(M)T~g(z,r p ), (2) 

where the factor of one half allows for the double counting of pairs, C 2S (At>,7) converts 
from redshift space to real space pairs, F gives the fraction of the pairs that will merge in 
the next T mg (the "last orbit" in-spiral time from r p ) and (M) is the mean incoming mass 
as estimated assuming a constant M/L. For the relatively massive galaxies considered here 
the dynamical friction is so strong that it is more violent relaxation with little timescale 
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dependence on the masses. The measured ratio of the numbers of 50 and lOO/i -1 kpc pairs 
to 20/i~ 1 kpc pairs in the CNOC2 sample is 3.8 ± 1.0 and 9.4 ± 3.0, respectively, in accord 
with the expectation of a growth as r^" 7 with the inner cutoff of 5/i _1 kpc. 

Not all kinematic pairs are close in physical space. The relation between the kinematic 
pairs closer than r p and Av and pairs with a 3D physical distance r p is readily evaluated by 
integrating the velocity convolved correlation function over velocity and projected radius 
and ratioing to the 3D integral of the correlation function. We find that C zs = 0.54 for 
Av = 1000 km s -1 and 7 = 1.8 There is support for this value on the basis of morphological 
classification, as tested in Patton et al. (1999), where about half of the kinematic pairs 
exhibited strong tidal features. 

The fraction of physically close pairs that are at sufficiently low velocity to merge is a key 
part of the rate calculation. It is clear that many galaxies will have close encounters which 
do not lead to immediate mergers, although mergers could of course occur on subsequent 
orbital passages. The key quantity that we need is the ratio of the critical velocity to 
merge, v mg , to o"i 2 . The timescale for close pairs to merge is much shorter than the time 
over which morphological disturbances are clearly evident, by nearly an order of magnitude 
(Barnes & Hernquist 1992; Mihos & Hernquist 1996; Dubinski Mihos & Hernquist 1999). 
This is one of our reasons for preferring kinematic pairs as a merger estimator. The 
simulation results indicate that the time to merge is, on the average, roughly that of a 
"half-circle" orbit, which at r p of 20/i _1 kpc at a velocity of 200 kms" 1 is close to 0.3 Gyr. 
A straight-line orbit with instantaneous merging would merge in about 0.1 Gyr, although 
that is not likely to be representative. 

To compute the merger probability, F(< v mg ), we need to know the maximal velocity 
to merge, v mg , at a physical separation of 20/i~ 1 kpc for a typical M* galaxy. A not very 
useful lower bound is fixed by the Keplerian escape velocity at 20/i _1 kpc, t> c v^2kms _1 , 
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where the circular velocity is approximately 200 kms -1 . An upper bound to v mg is the 
velocity that an object would have if it is captured into a galaxy's extended dark halo at the 
virialization radius and orbits to 20h~ 1 kpc with no dynamical friction. The virialization 
radius is approximately at the radius where the mean interior overdensity is 200p c , implying 
^200 = v c /(10H ), or, about 200/i~ 1 kpc for our typical galaxy. The largest possible 
apogalactic velocity at r 2 oo is v c , which leads to an undissipated velocity at 20h~ 1 kpc of 
2.37v c . Using a 12 = 200(300) km s' 1 at 20/i~ 1 kpc, we find that the fraction of all physical 
pairs that merge in one T mg is about 0.40(0.16). Therefore, we will normalize to a merger 
probability of 0.3, noting the 50% or so uncertainty. 

The absolute magnitude limit of —19.8 + 51og/i mag corresponds to L > 0.5L* which 
contains about 58% of the luminosity for the mean CNOC2 LF, M* = —20.4 and a = —1.2. 
To make our merger rate inclusive of major mergers we normalize to L > 0.2L*, which 
includes 85% of the luminosity. Within the current statistical accuracy, the paired and field 
galaxies have identical LFs. On the basis of n-body experiments (Barnes & Hernquist 1992) 
galaxies with masses greater than about 0.2M* will merge in approximately one orbital 
time. 

On the basis of these considerations we find that the rate of mass accumulation of 
galaxies with luminosities of 0.2M* and above is, 

TZ M = (0.02 ± 0.01)M*(1 + ^)0.i±o.5 ^m>i2) 0.3Gyr Gyr _^ (g) 

where we have adopted the 100/i~ 1 kpc mi value for a flat, low density cosmol- 
ogy and explicitly assumed that the velocity and timescale factors do not vary 
over this redshift range, as expected at these small scales in a low Q universe 
(Colin Carlberg & Couchman 1997). There is direct evidence that once evolution com- 
pensated that the luminous galaxies retain have no evolution in their circular velocities 
(Vogt et al. 1997; Mallen-Ornelas Lilly Crampton & Schade 1999). 
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5. Discussion and Conclusions 

Our main observational result is that for galaxies with M^' e < —19.8 + 5\ogh mag, 
the fraction of galaxy luminosity in 50/i _1 kpc wide kinematic pairs is about 14%, with 
no noticeable redshift dependence over the redshift zero to one range. This implies an 
integrated mass accretion rate of about 2% of per Gyr per galaxy for merging galaxies 
having L > 0.2L*. Our rate is uncertain at about the factor of two level due to uncertainty 
in the dynamical details of merging for our sample definitions. This merger rate implies a 
15% mass increase in an M* galaxy since redshift one. If the correlations of lower luminosity 
galaxies are only somewhat weaker than these (Carlberg et al. 1998) then the same 0.15M* 
merged- in mass causes a 50% mass increase in a 0.3M* galaxy. 

There are several issues that require further investigation. First, the rate of merging 
of similarly selected kinematic pairs should be studied in appropriately matched n-body 
experiments to better determine the orbital timescales. Second, the absence of a redshift 
dependence of a 12 and v mg needs to be observationally checked. Third, the connection 
between close kinematic pairs and morphologically disturbed galaxies, which does conform 
to the kinematic pair predictions at low redshift (Patton et al. 2000a), needs to be better 
understood at high redshift. 

This research was supported by NSERC and NRC of Canada. HL and DWH 
acknowledge support provided by NASA through Hubble Fellowship grants HF-01110.01- 
98A and HF-01093.01-97A, respectively, awarded by the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., for 
NASA under contract NAS 5-26555. 
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Fig. I. — The fraction of the sample luminosity in pairs with Av < 1000 km s -1 and 5 < 
r p < 50 hr 1 kpc as a function of redshift. The octagons are from the CNOC2 sample and 
the triangles from CFGRS. 

Fig. 2. — The velocity distribution function for CNOC2 galaxy pairs. The filled points are 
for 20/i~ 1 kpc pairs limited at M^' e of —19.5 mag (diamonds) and —18.5 mag (octagons) 
and the open squares are for 50/i _1 kpc pairs limited at —19.5 mag. The lines are the 
distributions expected for an r~ 18 correlation function convolved with Gaussian velocity 
distribution functions of 200 and 300 kms -1 , the upper and lower lines respectively. 
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